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Abstract 

Precise four-body calculations for ^He and have been performed in the 

framework of the variational method with Jacobian-coordinate Gaussian-basis 

functions. All the rearrangement channels of both NNNA and NNNT, are 

explicitly taken into account for the first time with the use of realistic A^A^ and 

YN interactions. The role of A-S conversion and the amount of the virtual 

S-component in ^He and are discussed. 
21.80.-ha,21.10.Dr,21.10.Gv,21.45.-Fv 
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I. INTRODUCTION 



It is one of the most fundamental problems in hypernuclear physics to extract novel 
information of YN interactions through precise calculations for a few-body systems such 
as ^H, ^He and ^H. Although a number of studies for ^He and have been performed 
with the use of various models so far, the relations to the underlying YN interactions 
are dependent on the adopted models In order to explore the features of the YN 

interactions, it is highly required to perform four-body calculation without any restriction 
in the configuration space. An especially interesting issue in this relation is to make clear 
the role of A-E conversion and to get a reliable estimate of the amount of the virtual S- 
component in A hypernuclei. 

However, the KN-HN coupling terms of the YN one-boson-exchange (OBE) models 
proposed so far have a lot of ambiguity due to scarce information from the YN scattering 
experiments. Recently, a series of the realistic YN interaction models called NSC97a~f ||5| 
have been proposed in which the spin-spin strengths are varied within the acceptable range in 
view of the limited YN scattering data. Correspondingly, by including the AN-'EN coupling 
explicitly, Miyagawa et al. performed the Faddeev calculations for extensively to test 
these OBE potential models in this lightest system. They found that the KN-HN coupling 
is crucial to get the bound state of and that only the version-f potential is acceptable 
among NSC97 models. For the further examination on the AA^-SA^ coupling, ^He and \R 
are much more useful because both of the spin-doublet states have been observed. Regarding 
the representative calculations done so far, Gibson et al. employed the coupled two-body 
model [|l| of '^IIe('^II)-|-A/S which was adopted originally by Dalitz and Downs 0, and 
then carried out the four-body coupled-channel calculation with the separable potentials 
of central nature 0. J. Carlson tried to perform four-body calclation of these hypernuclei 
with NSC89 model by using Variational Monte Carlo method ||^ and calculated the binding 
energies with statistical errors of 100 keV. Akaishi et al. 0] recently analyzed the role of 
the KN-TjN coupling for the 0^-1^ splitting also in the framework of the coupled two-body 
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model of ^He+A/S. Now our concern is to perform precise four-body calculations of NNNA 
and NNNH irrespectively of any model assumptions. 

It should be noted that, when one goes from to the A = A strange systems by 
allowing the A-S conversion, the computational difficulty increases tremendously. Recently 
we have successfully performed the extensive four-body calculations without any restriction 
on the configuration space: Both the NNNA and NNNT, channels have been incorporated 
explicitly and all the rearrangement channels of these baryons are taken into account. The 
variational method with the use of Jacobian-coordinate Gaussian-basis functions is 
adopted here; it has been proved to provide us with precise computational results for few- 
body systems |]8|-|T2|. 



The main purpose of this work is, first, to solve four-body problem of ^He and by 
taking into account the NNNA^E) channels explicitly with the use of realistic A^A^ and 
yA^ interactions and, secondly, to clarify the role of the AN-T,N coupling in the A = 4 
hypernuclei quantitatively. 

In Sec. II, we describe our method to solve the four-body problem and YN and A^A^ 
interactions adopted here. In Sec. Ill, we calculate the A binding energies, the E-mixing 
probabilities and the AA^ and SA^ correlation functions in ^He and \B.. The role of the 
AA^-SA^ coupling is investigated by dividing its whole contribution into the AA^ two-body 
and ANN three-body parts. A summary is given in Sec. IV. 

II. METHOD AND INTERACTIONS 

As the first step before going to the use of sophisticated OBE models, we employ here the 
AA^-SA^ coupled YN potential with central, spin-orbit and tensor terms [0 which simulates 



the scattering phase shifts given by NSC97f. The potentials parameters for central, spin- 
orbit and tensor terms are listed in Table I. 

The main reason of using this simulated version of NSC97f is to focus our attention 
clearly on the physical ingredients as well as for computational tractability. In this relation, 
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our most important criterion for selecting the YN interaction is that the observed binding 
energy of is reproduced reasonably: The interaction used here leads to the A-binding 
energy, i?A(^H), 0.19 MeV which agrees well with the observed data (i?A(|H) = 0.13 ± 0.05 
MeV). As for the A^A^ interaction, we employ the AV8 potential [jl4[. We shall, however, 
examine how the choice of the A^A^ potential affect the A binding energy; we additionally 
take the Bonn A potential [|l5l and the Minnesota potential (central force only) [|l^ and 



investigate the role of the tensor term of A^A^ potential. 

The total four-body wavefunction is described as a sum of the amplitudes of all the 
rearrangement channels (c = 1 — 4) of Fig. 1 in the LS coupling scheme: 

vl/,M(iHe,iH)= E EE E C^L'sw 

y=A,Sc=l al ss'Stt' 
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X A { [$ (r„ R„ pj [[Xs' (12)xi (3)].Xi {Y)], 

x[[r/,,(12)r/i(3)],r^,,(F)]^^J, (1) 
where the spatial part is expressed, with a set of quantum numbers a = {nl, NL, K, uX}, by 

$„/Af(r,R, p) = [0„i(r)^ArL(R)]K^j,A(p) 



IM 



Here, A is the three-nucleon antisymmetrization operator and x's and ?7's are the spin and 
isospin functions, respectively, with the isospin ty = (1) for Y = A (S). Functional form 
of 0nim(r) is taken as 0„irn(r) = r'- e~^^^^"^'^Yim{r) where the Gaussian range parameters are 
chosen to lie in a geometrical progression (r„ = ria"~^;n = 1 ~ Timax), and similarly for 
iPnlCR) and iu\{p)- These basis functions were verified to be suited for describing both 
the short-range correlations and the long-range tail behaviour of few-body systems 



Eigenenergies of the Hamiltonian and coefficients C's are determined by the Rayleigh-Ritz 
variational method. The angular momentum space of /, L, A < 2 is found to be enough to 
get sufficient convergence of the calculated results mentioned below. 
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III. RESULTS AND DISCUSSION 



All the calculations have been performed both for |He and ^H. Calculated of the 0+ 
ground state and the 1"*" excited state of |He and |H are illustrated in Fig. 2 together with 
the observed values (also see Table I). In the case of taking only the NNNA channel, both 
of the two states are unbound. Here, the NNNT, sector is divided into the (NNN)!!! and 
{NNN)3T, channels in which the three nucleons are coupled to isospin t — 1/2 and 3/2, 
respectively. When the (A^A^A^)iE channel is included, the 0"*" state becomes bound, but 
the 1+ state is still unbound. Then we found that the 1+ state becomes bound only when 
the {NNN)^T, channel is switched on. It is noted that the binding energy of the 0+ state 
increases only slightly with this t = 3/2 channel. The E-channel components turn out to 
play an essential role in the binding mechanism of the A — A hypernuclei, the [NNN)3ll 
channel being specially important in the 1+ state. The calculated binding energy of the 0"*" 
state almost reproduces the observed binding energy, while the 1+ state is less bound by 0.6 
(0.4) MeV for ^^H (iH), and hence the 0+-1+ sphtting is larger than the observed splitting. 
The calculated value of B^{\}le{Q+)) - 5A(iH(0+)) = -0.05 MeV is different from the 
experimental one, +0.35 MeV, although the Coulomb potentials between charged particles 
(p, 2=*=) are included. This difference should be attributed to the charge-symmetry-breaking 
component which is not included in our adopted YN interaction. 

Here, we examine how the calculated A binding energies depend on the properties of 
adopted NN potentials. Bando et al. once pointed out that the rearrangement effect 
originated from the A^A^ tensor force plays an important role for the mass dependence of 
A binding energies of yl = 3 ~ 5 hypernuclei. Now let us focus our concern on the A^A^ 
tensor forces: The three A" A" potentials compared here are the AV8 potential with a stronger 
tensor force, the Bonn A potential with a weaker tensor force and the Minnesota potential 
with a central force only. The binding energy of the |^H is reproduced reasonably when 
these A"A^ potentials are used together with our YN potential. In Fig. 3 the calculated A 
binding energies of the 0+ and 1+ states in ^He are shown in the cases of using these A" A" 
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potentials. It should be noted here that the larger ratio of the NN tensor component turns 
out to bring about the less binding of A. In the cases of using AV8 and Bonn A, in which the 
tensor compornents are included more or less realistically, the observed A binding energy 
is reproduced well. However, the use of the effective model with no tensor force such as 
the Minnesota potential gives rise to the over binding of A by 0.4 MeV. Hereafter, we report 
only the results for the AV8 potential, because there is no meaningful differece between the 
results for the AV8 and Bonn A potentials. 

As listed in Table HI, the calculated probabilities of the A'"A'"A'"E-channel admixture 
are 2.08% and 1.03% for the 0"*" and 1"*" states in ^He, respectively. In the 0"*" state, the 
probability of the (NNN)!!! channel is much larger than that of the (NNN)3'E channel, 
while in the I"*" state they are nearly the same. We therefore confirm that the {NNN)3Yi 
channel is especially important in the l"*" state. The S-, P- and D-state probabilities of 
the channels are also listed in Table HI. It is remarkable that, in the NNNT, channel, the 
£)-state component is dominant both in the 0+ and 1+ states, since the AN-T,N couphng 
part of the present interaction is dominated by the tensor component. These properties are 
almost similar in the case of ^H. 

It is of interest to see spatial locations of N, A and E particles in the A — A hypernu- 
clei and the correrations between each pairs. For the 0+ state of ^He, we illustrated the 
correlation functions (two-body densities) of the NN, AN and UN pairs (Fig. 4) and the 
(one-body) densities of single nucleon and A and E hyperons (Fig. 5). Also, we calculated 
the r.m.s. distances between the pair particles (fATTv? "^AAf, "^SAf) as well as the r.m.s. radii 
{"^N, ^A, ^s) of A^, A and E measured from the cm. of 3N (Table II). In Fig. 4, the NN corre- 
lation function in ^^He exhibits almost the same shape as that in the ^He nucleus, indicating 
that the dynamical change of the nuclear size due to the A participation is small. The AA^ 
correlation function is of larger range and flatter than the NN one, because the strength of 
the AN interaction is significantly smaller than the A^A^ case. The EA" correlation function 
is much shorter-ranged than the AA^ one due to the large virtual excitation energy (80 MeV) 
of A ^ E. These features are verified by the r.m.s. distances {fNN, ^aat, tsn) listed in Table 
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I. 

Furthermore, we notice that the A particle is located much outside the core nucleons 
(see fA > fiv and Fig. 5) and therefore the dynamical change of the core nucleus due 
to the A-particle partition is small: The nucleon r.m.s. radius is tn — 1-65 fm which is 
shrinked by 8 % from the corresponding one, — 1.79 fm in ^He. On the other hand, the 
E hyperon comes close to the nucleons (see tsat < taat and Fig. 4) and therefore generates 
a large dynamical contraction of the core nucleus in the A'"A'"A'"E-channel space; tn — 1-49 
fm is obtained with the E-channel amplitude only and the reduction rate amounts to 17 %. 
Furthermore, it is interesting to see in Fig. 4 the following feature of the E-admixture effect: 
In spite of the totally small probability of the E-mixing (2 %), the EA^ components at short 
distances are not so small in comparison with the AA" ones. This enhanced short-ranged 
component of the E-mixing is expected to be reflected in the non-mesic decay of EA^ — > A'A^. 

In order to show the physical effect of E-mixing in more detail, let us separate the whole 
contribution of the AN-T,N coupling interaction into the following two processes illustrated 
in Fig. 6: The first one is the process (i) which can be renormalized into the effective 
AA" two-body force and the second one is the process (ii) which can be represented by 
the effective ANN three-body force acting in the NNNA space. We solve the Schrodinger 
equation by excluding the three-body process (ii) so as to evaluate the contribution of the 
process (i) only. As shown in Fig. 6, the process (i) is large enough to make both the 0+ 
and l"*" states bound. We found that the three-body process(ii) is substantial: This process 
results in attraction by 0.62 MeV in the O"*" state, while repulsion by 0.09 MeV in the l"*" 
state of ^He and similarly for ^H. 

In order to investigate the respective role of the central and the tensor terms of the AA^- 
EA^ coupling interaction in the process (ii), let us make additional calculations with artificial 
modification as follows: We replace the AN-T^N coupling part of our original interaction 
by the purely central interaction and by the purely tensor one without changing the other 
interaction terms, where the calculated AA^ phase shifts are kept almost unchanged. It is 
found that the effect of the process(ii) is quite different between the two cases. Namely, in 
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the case of the purely central-force coupling, the contribution is —1.54 MeV (attractive) in 
the 0"*" state and +0.43 MeV (repulsive) in the 1"^ state. On the other hand, in the case 
of the purely tensor-force coupling, it becomes much reduced values, —0.35 MeV for the O"*" 
state and —0.06 MeV for the I"*" states. Thus, the resulting energies of 0"^ and 1+ states 
turn out to be controlled substantially by the ratio of the central and the tensor components 
in the AN-T,N coupling potential. Characteristic feature of our result of the process (ii) in 
the 0"*" state is similar to that given by Akaishi a/ @]. In the 1"*" state, however, the effect 
of the process (ii) is much different from the result of Ref . [Q : In our treatment the process 
(ii) works repulsively or attractively depending on the ratio of the central and the tensor 
components in the AA^-SA^ coupling potential, while the contribution of the process (ii) is 
negligible in Ref. 0] independently of it. This difference may come from the fact that in the 
1^ state of aHc, the probability of the {NNN)3'E is as large as that of the {NNN)iT, in our 
result, but the (NNN)!!] channel cannot be explicitly included in the coupled two-body 
model of =^He + A/S of Ref. g. 

IV. SUMMARY 

We have developed the calculational method of four-body bound-state problems so that 
it becomes possible to make precise four-body calculations of and ^He taking both the 
NNNA and NNNT, channels explicitly into account using realistic A^A^ and YN interac- 
tions. As a result, we succeeded in making clear the role of A-S conversion and deriving 
the amount of the S-mixing in A = 4 hypernuclei quantitatively. However, our YN inter- 
action employed here is not sufficient to reproduce the binding energy of the excited state 
of 1"^, although those of ground states of |H, ^He and are in good agreement with the 
experimental values. It is a future problem to explore the feature of A-S conversion in 
A hypernuclei with the use of more refined YN interactions through systematic study of 
structure of the heavier hypernuclear systems. 
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Table I. Parameters of the central (C), spin-orbit (LS) and tensor (T) terms of the 
NSC97f-simulated YN interaction. Radial form of each term is given by two-range Gaus- 
sians. Range is in units of fm and strength (V) in units of MeV. 





P 






Vls 


Vt 


NK-NK 


0.5 


732.08 


1068.8 


1023.8 


-243.31 




1.2 


-99.494 


-45.490 


-17.195 


-10.413 


NK-NT, 


0.5 


61.223 


-770.21 


-19.930 


287.54 




1.2 


-15.977 


68.274 


22.299 


62.438 


nj:-nj: 


0.5 


1708.0 


863.76 


544.56 


21.778 


(t = 1/2) 


1.2 


80.763 


28.284 


-19.944 


-53.542 




0.5 


695.39 


-181.08 


-462.31 


333.05 


{t = 3/2) 


1.2 


-109.37 


23.282 


0.0023 


22.234 
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Table II. Calculated energies of the 0+ and 1+ states of ^^He and ^^H. The energies E are 
measured from the NNNA four-body breakup threshold, va-b denotes the r.m.s. distances 
between particles A and B, while ta stands for the r.m.s. radius of particle A measured from 
the cm. of 3A^. As for the ^He (^H) nucleus, the calculated binding energy is —7.12 (—7.77) 
MeV and fN-N = 3.10 (3.03) fm and fjv = 1.79 (1.75) fm. 



J 


0+ 


1+ 


0+ 


1+ 


E{MeV) 


-9.40 


-7.66 


-10.10 


-8.36 


E'^^P(MeV) 


-10.11 


-8.87 


-10.52 


-9.53 


5a (MeV) 


2.28 


0.54 


2.33 


0.59 


Sl"P(MeV) 


2.39 


1.15 


2.04 


1.05 


fN-N (fm) 


2.86 


3.03 


2.83 


2.99 


fA-N (fm) 


3.77 


5.74 


3.75 


5.70 


fs_Ar (fm) 


2.24 


2.48 


2.23 


2.46 


Tn (fm) 


1.65 


1.75 


1.64 


1.73 


fA (fm) 


3.39 


5.47 


3.37 


5.43 


fs (fm) 


1.67 


1.81 


1.66 


1.80 
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Table III. The S-, P- and D-state and total probabilities of the {NNN)iA, {NNN)iJ: 
and {NNN)3J: channels in the 0+ and 1+ states of (a)^He and ^H. {NNN)t denotes three 
nucleons whose isospins are coupled to t. 

(a) iHe 



(b) iH 





s (%) 


P(%) 


D{%) 


Total (%) 


0+ 

{NNN)iA 


89.32 


0.08 


8.52 


97.92 


{NNN)iJ: 


0.84 


0.04 


1.16 


2.04 


{nnn)3j: 


0.01 


0.01 


0.02 


0.04 


1+ 

{NNN)iA 


90.38 


0.07 


8.52 


98.97 


{nnn)ij: 


0.10 


0.01 


0.40 


0.51 


{nnn)3j: 


0.09 


0.00 


0.43 


0.52 






S (%) 


P (%) 


D (%) 


Total (%) 


0+ 

{NNN)iA 


89.27 


0.08 


8.54 


97.88 


{NNN)iJ: 


0.86 


0.04 


1.19 


2.08 


{nnn)^j: 


0.01 


0.01 


0.02 


0.04 


1+ 

{NNN)iA 


90.34 


0.07 


8.54 


98.95 


{nnn)ij: 


0.10 


0.01 


0.41 


0.52 


{nnn)3j: 


0.09 


0.00 


0.44 


0.53 
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Figure caption 

Fig. 1 Jacobian coordinates for all the rearrangement channels of NNNAiT,) system. 
Three nucleons are to be antisymmetrized. 

Fig. 2 Calculated energy levels of (a)^He and (b)^H. The channels successively included 
are (i) {NNN)iA, (ii) (iViViV) i E and (iii) (iViViV) 3 E where the isospin of the three nucleons 
is coupled to t =1/2 or 3/2. Energy is measured from the ^He+A(^H+A) threshold. 

Fig. 3 Calculated energy levels of ^He with the AV8, Bonn A and Minnesota potentilas. 

Fig. 4 Correlation functions (two-body densities) of the NN, AN and EA^ pairs in the 
O"*" state of ^He together with that for the NN pair in ^He. Here, the correlation function of 
T,N pair has been multiphed by factor 2 so as to see the behavior of this function clearly. 

Fig. 5 Calculated one-body ddensities of A^, A and E particles in the 0+ state off ^He. 
Volume integrals of the densities are 1.0, 0.98 and 0.02 for A^, A and E particles, respectively. 

Fig. 6 Calculated energy levels of ^He for case (i) and case (i)-|-(ii); here, (i) denotes 
the two-body process and (ii) denotes the three-body process. 
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